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We propose the model of first quasars formation around the cluster of primordial black holes (PBHs). It is supposed, 
that mass fraction of the universe ~ 10~^ is composed of the compact clusters of PBHs, produced during the phase 
transitions in the early universe. The clusters of PBHs become the centers of dark matter condensation. As a result, 
the galaxies with massive central black holes are formed. In the process of galaxies formation the central black holes 
are growing due to accretion. This accretion is acompaned by the early quasar activity. 
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1. Introduction 

The discovery of distant quasars with the redshifts 
2; > 6 in the Sloan Digital Sky Survey P provides new 
questions in the problem of galaxy formation. The 
maximum observed quasar redshift z = 6.41 corre- 
sponds to corresponds to accretion on the black hole 
(BH) with the mass S-IO^Mq gj. The early formation 
of BHs with masses ~ lO^M© puts serious difficulties 
on the BH formation scenario due to the dissipative 
evolution of star clusters 0], supermassive or gaseous 
disks 0] . In addition it is dfficult to reconcile the shape 
of the quasar redshift distribution with the hypothe- 
ses of slow continuous BH growth in the processes of 
their merging or accretion j5| 1^ . For these reasons the 
scenario of PBHs origin j7| IS] becomes more attrac- 
tive. These PBHs can be the centers of baryon "5 and 
dark matter (DM) JO] condensations into the growing 
galaxies. 

The new effective mechanism of PBHs formation 
was developed in ^2 ^1 . This model predicts the 
specific cluster structure and properties of the forming 
PBHs ^3], In the current paper we will use this 
model for describing the formation of galaxies around 
the clusters of PBHs. It is the clusters of PBH who 
could play the role of the initial density perturbations. 
As the basic example, a scalar field with the tilted Mex- 
ican hat potential was accepted. Properties of PBHs 
clusters appear to be strongly dependent on the ini- 
tial phase of the scalar field. In addition they strongly 
depend on the tilt of the potential and the scale of 
symmetry breaking / at the beginning of the inflation. 
Here we elaborate the dynamics of DM coupled with a 
PBH cluster by gravity. It is shown that the galaxies 
could be formed in this case even in the absence of the 
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Figure 1: The initial mass profile Mh{ri) of PBHs 
cluster and the DM mass profile AfDM(f'i) • 



primordial fluctuations in the DM. We will use here 
the same values of parameters as in the j^lE], which 
are in the agreement with observations. 

The initial mass proflle Mh{ri) of PBHs cluster is 
calculated here according to 14 and has the speciflc 
form presented in the Fig. ^ For comparison the mass 
of DM M-£,M{ri) inside the same spherical shells are 
shown. Radius denotes the physical size at the mo- 
ment ti and the temperature Ti , when the correspond- 
ing sphere are crossing the cosmological horizon. Any 
physical size at the temperature Ti is smaller than that 
in the recent epoch at T^/Ti times, where Tq = 2.7 K. 
Note that shells in the Fig.^are taken at different mo- 
ments ti . Respectively the shown mass of uniformly 
distributed DM does not follow to the law A/dm oc 
as it must be for flxed time case. 

The mass distribution of PBHs is strongly distorted 
after their formation due to the fractal structure of 
the cluster. Indeed, an each PBH in the cluster is ac- 
companied by the subcluster of smaller mass PBHs. 
The much more massive PBH is formed after the cor- 
responding gravitational radius radius crosses the hori- 
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zon. The density in the local center is so high that a 
lot of PBHs appear to be inside their common gravita- 
tional radius Vg = 2GM /c^ . The total mass appears 
inside the horizon is 4.3- lO^Af© . So the whole range of 
masses and radii shown in the Fig.^is not realized. As 
a result we obtain the initial mass distribution of PBHs 
with the much more massive PBHs than it could be ex- 
pected. In fact, the satellite PBHs of smaller masses 
increase this value in several orders of magnitude. The 
next Sections are devoted to the investigation of mu- 
tual evolution of DM (uniformly distributed from the 
beginning) and the cluster of PBHs with distribution 
according to the Fig. ^ Note that the PBHs them- 
selves contribute to the DM, but we don't consider 
them as the main DM part. 

There are several stages of PBHs and galaxies for- 
mation: (i) The formation of the closed walls of scalar 
field just after inflation and their collapse into the clus- 
ter of PBHs according to . Formation of the most 
massive BH in the center of the cluster after the hori- 
zon crossing, (ii) The detaching of the central dense re- 
gion of the PBHs cluster from the cosmological expan- 
sion and its virialization. Many of surrounding small 
BHs merge with a central most massive BH and in- 
crease its mass. (iii). The quasars ignition due to ac- 
cretion onto the central BH. (iv) Detaching of the outer 
cluster region where DM dominates from the cosmolog- 
ical expansion and further galaxy growth. Stop of the 
galaxy growth due to interaction with the surrounding 
DM fluctuations originated from the inflation, (v) The 
gas cooling and star formation. 

There are definite astrophysical limitations on the 
number and mass of PBHs: The universe age limit 
gives for PBH cosmological density parameter fi/i < 1 . 
The possibility of tidal destruction of globular clusters 
by the PBHs gives the PBH mass Hmit Mh < IO^Mq 
if these PBHs provide the major contribution to the 
DM . The contribution of accreted PBHs into the 
background radiations gives approximately flh < 10^^ 
for Mh ~ IO^Mq [m. The limit Qh < 0.01 for 
lO^M© < Mh < lO^Mo is obtained from the VLBI 
observations of the lensing of compact radio-sources 
|17j . In this paper we consider the case flh ^ 10^^ 
which is in accordance with all aforementioned limits. 

In the following the subscript '0' marks the values 
of a current time tg , 'eq' corresponds to the time t^q of 
matter-radiation equality, and 'i' to the time of horizon 
crossing respectively. We consider the flat cosmological 
model with the density parameters ^m.o = 0.3, f^A.o = 
0.66, Qb.oh'^ = 0.02 and h = 0.7. 



2. Gravitational dynamics of cluster 

Let us consider the gravitational dynamics of PBHs 
cluster and the DM. The results of this section are ap- 
plicable as for the inner part of the cluster, composed 



mainly of PBHs and collapsing at radiation dominated 
stage, and for outer regions of the cluster, where the 
DM is the main dynamical component. The later re- 
gions are detached from the cosmological expansion at 
the matter dominated epoch. The spherical symme- 
try is supposed. This approximation is rather good for 
the inner regions because the considered density fluc- 
tuations have a large amplitude in comparison with the 
standard inflationary generated fluctuations. 

Consider a spherically symmetric system with the 
radius r < ct, consisting of PBHs with the total mass 
Mh inside the radius r , the radiation density pr , the 
DM density pdm and the cosmological A -term density 
Pa . The radiation density (and obviously the density 
of A-term) is homogenous. Therefore, the fluctuations 
induced by the PBHs are of the type of entropy fluc- 
tuations. Because the scale under consideration is less 
then the cosmological horizon scale we use the Newto- 
nian gravity but take into account the prescription of 
jl8j to treat the gravitation of homogenous relativistic 
components p p + 3p/c^ . the evolution of spherical 
shells obeys the following equation 



G{Mh + Mdm) SirGprT SwGpAr 



(1) 



with the approximate initial conditions at the moment 
U: r — —Hr, r{ti) = . During the derivation of 
Eq. ^ it was taken into account, that Sr + 3pr — 
2er, £a + SpA = — 2eA- In the parametrization r = 
£,a{t)b{t) , the ^ is the comoving length, a{t) is the 
scale factor of the universe normalized to the present 
moment as a(io) — 1 and the function b{t) char- 
acterizes the deflection of the cosmological expansion 
from Hubble law. The ^ connected to the mass of 
DM inside considered spherical volume (excluding BHs 
mass) by the relation Mdm — (47r/3)pDM(to)'?'^ • Func- 
tion a{t) obeys the Friedman equation 



SttG . , 

— 5-(£r + £m + Ea), 

3c- 



(2) 



which can be rewritten as aja = H()E{z), where red- 



shift z 



- 1 , Hq is the current value of the Hubble 



constant and function 

E{z) = [QrA^ + zf + r!„,,o(l + zf + ^AflY'^. (3) 

Here VLr , and J^a is the density parameters (frac- 
tions of density to critical density pc = 3H'^ / {8ttG)) 
of radiation (in sum with relativistic neutrino), non- 
relativistic matter (DM and baryons) and A-term, re- 
spectively. The A-term affects the evolution of per- 
turbations only at the redshift z < 10. Relations t{z) 
and z{t) can be easily obtained from the solution of 
the equation 

The considered cosmological model is flat: ilr + 
r2„i + CIa = 1 • By using the second Friedman equation 
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Figure 2: The virial radius of galaxy is shown as a 
function of redshift z . 



(for a) one can rewrite as follows 



db ^ 
dz 



52 



-b 



^m,o(l + z) 

2E^{z) 



where function 

1 dE{z) 



S{z) = 



E{z) dz 



1 



0, (4) 



(5) 



and value of fluctuation 5h — Mh/M-QM- In the limit- 
ing case ea = the is equivalent with the equation, 
which was obtained at JHI- Unfortunately, the fitting 
formula of ^ p ~ 140$^($ + l)poq, where <^ = 5h, 
does not allow us to consider the very dense clusters 
(with (5/j > 10'') and we must solve Eq. 10}. In de- 
pendence of epoch one can omit some terms in the 
r.h.s. of iQJ. We start tracing the evolution of cluster 
at high redshift Zi when the considered shell crosses 
the horizon r ^ ct. Initial conditions for evolution are 
presented in the Fig. ^ 

The moment of the expansion stop f = is defined 
by the condition db/dz = b/{l + z). We put that 
after expansion stop. After contraction from maximum 
radius to rc = rs/2 the shell is virialized and fixed 
at the radius Tc . Therefore, the mean density of the 
DM in the virialized object p is 8 times greater then 
the DM density at the moment of maximum expansion 



P = 8p™,o(l + Zsfbg ^ 
and an effective (virial) radius of the object 



/3Mdm\ 
V 47rp ) 



1/3 



(6) 



(7) 



The radius Q is shown in the Fig. [2 

Let us consider the fate of spherical shells by mov- 
ing our attention from the center of the cluster to outer 
regions. It is obvious (and supported by numerical cal- 
culations) that inner more dense shells stop expansion 
early than outer shells. As was discussed before, a BH 
with the mass Mc = 4.3 • lO^M© forms in the center of 
the cluster at the moment ti . The subsequent shells. 
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Figure 3: The final density profile © of a galaxy (p 
in units Mq/pc^) in dependence of the distance from 
the center of cluster rc is shown for DM (dashed line), 
for BHs (dotted line) and for the sum of these densities 
(solid line). Asymptotic power laws are also marked. 



where > 1 (Mdm < Af/i), are very dense too. These 
shells are detached from the cosmological expansion at 
the radiation dominated epoch. In general, these shells 
accreted by the central BH in the process of two-body 
relaxation of BHs. This effect will be considered be- 
low. The boundary value 6h — ^ corresponds to the 
masses Mh = Mdm = 3.3 • 10*Mq . 

For the outer shells where the fluctuation is small 
(5/,. < 1 , the known Meszaros solution is true, according 
to which the fluctuation growth till the moment teq is 
equal to 2.5. For the early formed PBHs it is possi- 
ble the process analogous to the 'secondary accretion'. 
As a result the PBHs would be 'enveloped' by the DM 
halo. We call these haloes as induced galaxies (IG). 
The density profile does not follow the secondary ac- 
cretion law p (X r~^/* because the central mass is not 
compact. After virialization the distribution of DM is 



PTiM{r) 



1 dM^M{r) 
Anr^ dr 



(8) 



where function Mdm(''c) is determined by the solution 
of 0). In analogy with the DM, one can obtain the pro- 
file of the BHs density and of total density. The corre- 
sponding results are shown in the Fig.|31 where density 
is expressed in units Mq/pc^ and distance from the 
center is in pes. As can be viewed from the Fig. |3| at 
the Sun distance from the Galaxy center r = 8 kpc the 
total local mass density is 0.7 Gev cm~^. Therefore, 
some characteristics of the obtained galaxy are similar 
to the our Galaxy. But the considered induced galaxy 
is more dense at its center and hosting the supermas- 
sive BH. 

The resultant structure of induced galaxy is the 
following. Inside radius r = 10 pc from the induced 
galaxy center the mass of BHs (including a central su- 
permassive BH) is 2.9 • IO^Mq . Inside the same radius 
the mass of DM is 3.4 • lO^M© . At the larger distances 
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the density profile of DM is well approximated by the 
power low 



p{r) = 2.2 • 10'' 



10 pc 



-2.2 



Mqpc-3. 



(9) 



The density profile © differs from the Navarro-Frenk- 
White [50] and other proposed profiles, obtained from 
the numerical simulation of haloes formation (without 
primordial cluster of BHs) but corresponds to these 
profiles at intermediate scales where power index « 
—2. An interesting properties is the diminishing of 
mean velocity K = {GM/2RY/^ of IGs in time (with 
decreasing of z). Such a behavior is a consequence 
of the shape of perturbation spectrum produced by 
clusters of BHs (see next Section). 

3. Termination of induced galaxies 
growth 

The growth of a virialized region terminates at the 
epoch of nonlinear growth of the ambient density fluc- 
tuations of the same mass M ( originated from the 
standard inflation ) as the our combined system of 
PBHs plus the DM halo. The laws of growth for the 
both fluctuations at the matter dominated epoch are 
the same (|13|) . Therefore, the condition of the growth 
termination of a typical induced galaxy is the equal- 
ity of r.m.s. perturbations (5gq(AfDM)7 produced by 
the PBHs cluster, and the inflationary r.m.s. pertur- 
bations (j^^{My)m) , both depending on the DM mass 
scale A'/dm: 



.DM 



(10) 



where v is the perturbations peak height, in this sec- 
tion we consider only mean perturbation with v — I . 
The statistical properties of the inflationary Gaussian 
perturbations are determined by the power spectrum 
P{k) . We use the power spectrum of the DM from 
EH: 



P{k) = 



Ak 



(1 + 1.71u + 9ui-5 + 1*2)2^ 



(11) 



where u = k /[{^m,o+^b,o)h'^ Mpc~"^] , and k is the co- 
moving wave vector in Mpc~^ units. The initial spec- 
trum is supposed to be the Harrison-Zeldovich spec- 
trum. The normalizing constant A corresponds to the 
value 0.9 of r.m.s. fluctuations at 8 Mpc scale. The 
r.m.s. perturbation at the mass scale M correspond- 
ing to the radius R is 



a^^{M{R)) = I k^dkP{k)W{k,R), (12) 



where W{k, R) is a filtering function. The evolution 
of the cosmological perturbations in the universe with 



PBHs cluster 



Inflation 
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Figure 4: The r.m.s. density perturbation values at 
the time tcq of matter-radiation equality produced by 
the cluster of PBHs and by the inflation respectively. 



the A-term at the matter dominated epoch can be 
obtained from or simply from the equation El] 

S{t)/S{zcq) = .9(2) (1 + Zcq)/[g(ioq)(l + z)] , whcrc 

g{z) « ^r!„jr!i/^-f}A+(l+r!„,/2)(l+r!A/70)]-i,(13) 



where il„i = ^^m,o(l + z)^ /E'^{z) and so on. 

The value of fluctuations produced by the cluster 
S'^^{Mum) = 2.5(5i(AfDM) is shown in the Fig. g] to- 
gether with the Harrison-Zeldovich spectrum (|ll|l in- 
serted to (|12f) . Fluctuations produced by the cluster 
are very well approximated by the power law 



(S^\(M) = 0.044 



M 



lOiOMp 



-0.85 



(14) 



From the numerical solution of (|1U|I we obtain the 
mass of the DM halo at the moment of growth termi- 
nation A/dm = 1-8 • IO-'^Mq . The expansion stops at 
the redshift z = 1.64. Note that growth termination 
does not mean the termination of large scale structure 
formation. This means only the termination of the 
contraction under the influence of PBHs cluster. After 
this moment the structure formation proceeds by the 
standard hierarchical clustering scenario: galaxies are 
assembled into larger ones, clusters and superclusters. 
Therefore, the large scale structures in the our model 
is the same as in the standard cosmological scenario. 

The formed induced galaxies looks like the giant 
elliptical galaxies with a small ellipticity. They have 
the central DM spike shown in the Fig. [3] and central 
supermassive BH. 

4. Evolution of quasar activity 

In view of the discovery of distant quasars with the red- 
shifts z > 6 , let us discuss the properties of induced 
galaxies at these redshifts. The described model of in- 
duced galaxies can easily explain the AIq — 3 • IO^Mq 
supermassive BHs at z — 6 (the age of the universe 
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= 9 • 10* year). In this model there are supermassive 
BHs formed at the radiation dominated epoch long be- 
fore the galaxies formation. Initially there is at least 
Mc — 4.3 • l(f Mq or even larger (if there is an adsorb- 
tion of surrounding smaller BHs from cluster by the 
most massive one) BH at the cluster center. Let us 
suppose that this BH radiated at the Eddington limit, 
Le = 1.3 • lO^^erg/s. This BH grows exponentially 
with a characteristic time Ie ^ 4.5 • 10*7^ yrs, where 
T] is the effectiveness of the accretion matter-energy 
transformation. The time of growth from Mc = 4.3 • 
IO^Mq to A/q ^ 3 • IO^Mq is Ai Ie In (Mq/M^) ~ 
2 • 10* year for ry ~ 0.1 . This means that QSO activity 
began at the moment ~ At ~ 7 • 10* yr, correspond- 
ing to the redshift z — 7.3 . At this redshift the induced 
galaxies have the virial mass 4.3 • IO^Mq , the virial 
radius 40kpc, the mean virial velocity 140 km/s and 
the mean virial temperature = rUpV^ /3 ~ 2 • 10^ K, 
where nip is the proton mass. These values were ob- 
tained from the numerical solution of the equation Q). 
For these parameters the radiative cooling mechanism 
for the gas in the induced galaxies is effective. This 
provide the cool gas flows necessary for the accretion. 

The very popular model of quasars ignitions is the 
gas flow onto the central black hole due to the tidal in- 
teractions during the merging of galaxies. In our model 
the huge tidal interactions are produced at the time of 
termination of the induced galaxies growth considered 
in the preceeding Section. In general, the density dis- 
tribution of surrounding inflationary perturbations are 
not spherically symmetric themselves and with respect 
to the induced galaxy center. This causes the strong 
tidal forces, instabilities and gas supply to the central 
region necessary for the effective accretion and quasar 
activity. We expect the maximum activity at the time 
of the termination of the induced galaxies growth at 
z ~ 1.6. This epoch is in agreement with the observ- 
able shape of quasar redshift distribution. 

For the more detailed description of induced galax- 
ies formation it is important to take into account dis- 
tributions both the induced galaxies (or initial PBHs 
clusters) and inflation perturbations. For simplicity, 
let us consider only the Gaussian distribution of infla- 
tion perturbations over v . In a, mean = 1 , the fluc- 
tuations produced by BHs cluster till z = 1.6 prevail 
over the standard inflation perturbations. But there 
are the less common fluctuations with v > I which 
exist at the discussed redshifts z — 6, 7.3 and larger. 
These rare fluctuations responsible in the our model for 
the tidal forces and quasar ignition at large z . From 
the numerical solution of I^J we can obtain the virial 
mass of induced galaxies in dependence of the redshift 
M-DMiz) ■ Substituting this function into H10() we get 
the dependency over ^{z) . The i^(z) is the required 
peak height of the inflationary fluctuation for the ter- 
mination of the induced galaxies growth and quasar ig- 
nition. The distribution of probabilities of these events 
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Figure 5: The form of quasar ignition events distribu- 
tion beyond the redshift z = 1.6 

is calculated as 

oo 

/(,) = _^ / 1 e-^'V^d.' = 1 e-'^W/^^ 
^ ^ dz J V2^ V2^ dz 

u(z) 

and shown in the Fig. [S] This Figure shows the dis- 
tribution law of the quasar ignition events beyond 
z = 1.6 where the hierarchical clustering begins. 
Note, that the Fig. |31 doesn't present quasar dis- 
tribution function Fqso{z). To calculate Fqso{z) 
we need the poorly deflned quantity — the time of 
quasar activity tq{z,MBH) (see for details). For a 
given tq{z, Mbh) it is possible to calculate -Fqso(^) — 
6z{tq, z)f{z) , where Sz(tq,z) is the redshift interval 
corresponding tq at redshift z . 

5. Conclusion 

In this paper we describe the new model of quasar for- 
mation initiated by the earlier formation of the clus- 
ter of PBHs. This model provides the early formation 
of large galaxies with the massive BHs in their cen- 
ters. Nowadays these galaxies could be seen as distant 
quasars. The main calculated parameters of typical 
galaxy are the following: the galaxy mass 2 • IO^^Mq , 
the central BH mass 4 • lO'^ - 10^ Mq . The density 
proflle of these induced galaxies is a near isothermal, 
p cx r~^-^ . The induced galaxies could be the cause of 
the of early quasar activity. The model predicts the 
rapid decay of the quasar activity at redshifts z > 10. 
The alternative scenario could be based on the forma- 
tion of the less massive PBHs clusters. In this case 
each nowadays galaxy contains numerous BHs origi- 
nated from the dwarf induced galaxies aggregated in a 
single galaxy in the process of hierarchical clustering. 
This possibility will be considered in a separate paper. 
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